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26Abstract 


27 Understanding how phytoplankton interacts with local and regional drivers as 
28well as their feedbacks is a great challenge, and quantitative analyses of the regulating 
29role of human activities and climate changes on these feedback loops are also limited. 
30By using monthly monitoring dataset (2000-2017) from Lake Taihu and empirical 
31dynamic modelling to construct causal networks, we quantified the strengths of causal 
32feedbacks among phytoplankton, local environments, zooplankton, meteorology as 
33well as global climate oscillation. Prevalent bidirectional causal linkages between 
34phytoplankton and the tested drivers were found, providing holistic and quantitative 
35evidence of the ubiquitous feedback loops. Phytoplankton exhibited the highest 
36feedbacks with total inorganic nitrogen and ammonia and the lowest with nitrate. The 
37feedbacks between phytoplankton and environmental factors from 2000 to 2017 could 
38be classified by two groups: the local environments (e.g., nutrients, pH, transparency, 
39zooplankton)-driven enhancement loops promoting the response of the phytoplankton, 
40and the climate (e.g., wind speed)-driven regulatory loops suppressing it. The two 
4lcounterbalance groups modified the emergent macroecological patterns. Our findings 
42revealed that the causal feedback networks loosened significantly after 2007 
43following nutrient loading reduction and unsuccessful biomanipulation restoration 
44attempts by stocking carp. The strength of enhancement loops underwent marked 
45decreases leading to reduced phytoplankton responses to the tested drivers, while the 
46climate (decreasing wind speed, warming winter)-driven regulatory loops increased- 


A7like a tug-of-war. To counteract the self-amplifying feedback loops, the present 


48eutrophication mitigation efforts, especially nutrient reduction, should be continued, 
49and introduction of alternative measures to indirectly regulate the critical components 


50(e.g., pH, Secchi depth, zooplankton biomass) of the loops would be beneficial. 


51 
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54Introduction 


55 Lake eutrophication, which cause a shift from clear to turbid water in shallow 
S6lakes , has become a long-standing environmental challenge to water quality 
57managers worldwide . This shift is largely attributed to an increase in external nutrient 
58loading (e.g., nitrogen, N and phosphorus, P), breaking the macrophyte-dominated 
59positive feedback loops in the clear water state and establishing phytoplankton 
60dominated positive feedbacks . In the clear water state, macrophytes prevent sediment 
6lresuspension, oxidise the sediment, enhance water transparency, outcompete 
62phytoplankton and eventually promote macrophyte growth , which lengthens and 
63suppresses internal nutrient recycling . Phytoplankton-dominated feedbacks might 
64form enhancement loops (e.g. elevated water pH and organic aggregates, reduced 
6S5dissolved oxygen near the sediment and low water transparency, and stimulated P 
66release from the sediment) that may lead to super linear relationships between 
67phytoplankton and nutrients , while macrophyte-associated feedbacks working as 
68important regulatory reducing phytoplankton proliferation (Moss, 1990; Søndergaard 
69et al, 2017). Thus, there is growing evidence that the enhancement feedback loop 
70would maintain and prolong the algae bloom ), which, in turn, may delay the lake 


71rehabilitation processes and counteract restoration efforts taken by the governments . 


72 As any variations in the components (i.e., phytoplankton, nutrients, pH, SD, 
73zooplankton) of the feedback loop might act as potential drivers enlarging or 
7Alessening its ecological consequences , external nutrient reduction is usually a 


75practical and feasible way to combat eutrophication as well as nuisance blooms . 


76Likewise, the feedback loops may also be depressed by alternative eutrophication 
77mitigation measures such as application of flocculants to decrease suspended 
78particular and particular P as well as chemical inactivation agents to preserve P in 
79sediments, and biomanipulation to control phytoplankton through a trophic cascade . 
80However, as a coin has two sides, a series of potential drivers may push the feedback 
8lloop into a worse state by altering its components . For example, climate warming 
82facilitates algae growth and aggravates the nuisance blooms , extreme rainfall brings 
83nutrients (N, P) and suspended particulars (decreased SD) into lakes , and atmospheric 
84stilling decreases non-algae turbidity and thus enhances the light availability for algae 
85production . Therefore, to revert the feedback loops, lake managers may have to spend 


86substantial investment for a long time period to combat deterioration . 


87 Take Lake Taihu as an example. It is the third largest lake in China, and 
88infamous for the water crisis in 2007 (Guo, 2007; Qin et al., 2010). Though a series of 
89intensive restoration efforts have progressively been implemented since 2007 (Qin et 
90al., 2019), the frequency and magnitude of algae blooms were not decreased as 
9lexcepted (Qin et al., 2021; Qin et al., 2019). For instance, in 2017, chlorophyll a (Chl) 
92increased by nearly 80% but total phosphorus (TP) only by 31% compared to the 
93average value across years (2005-2019) , suggesting super-linear (s/ope >1) relations 
94(macroecological patterns) of logChl-logTP. This phenomenon might be attributed to 
95a dominant self-amplifying feedback loops that regenerates nutrients to fuel algae. 
96While the intensive restoration efforts might diminish the enhancement feedback 


97loops or strengthen the regulatory feedback loops which could relieve cyanobacterial 
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98bloom (Qin et al., 2019; Schindler et al., 2016; Smith and Schindler, 2009), climate 
99changes (e.g., warming, extreme rainfall, atmospheric stilling) might impede it (Yang 
100et al., 2016; Zou et al., 2022), resembling a tug-of war (Guo et al., 2019; Meerhoff et 
10lal., 2022; Qin et al., 2021; Yan et al., 2017). So far, a holistic and quantitative 
102assessment of the existence and strength of these self-reinforcing feedback loops as 


103well as how they are modified by external drivers is unclear. 


104 In this study, we assembled a long-term (2000-2017) dataset on Lake Taihu to 
105construct the causal networks involving how phytoplankton interacted with water 
106nutrients (i.e., nitrate, NO3, ammonium, NH4, total inorganic nitrogen, TIN and 
107orthophosphate, PO4), pH, SD and zooplankton. With the dataset we tested for the 
108following hypothesised feedback loops (Fig. 1): nutrients > Chl (1.e., chlorophyll a, 
109used as a proxy of phytoplankton biomass) <> zooplankton <> nutrients; nutrients +> 
110Chl <> pH/SD < nutrients and referred it to the emerging macroecological patterns 
111(EMPs, the estimated slope of the log-log linear relationships) of the phytoplankton 
112response to bottom-up (e.g., TP, TN) and top-down (e.g., zooplankton) controls. We 
113were also interested in exploring the causal linkages between the components of the 
114local feedback loops and the regional and global drivers. Specifically, we addressed 
115the following questions: (1) What is the relative importance of the distinct causal 
116determinant for Chl? (2) which are the predominant local feedback loops associated 
117with Chl and what is the potential influence of global climate oscillation and the 
118regional meteorological forces on the components of feedback loops within the lake? 


119(3) what is the relative importance of the enhancement and regulatory loops for 


120determining the EMPs of the phytoplankton response to bottom-up (e.g., TP, TN, SD) 
12land top-down (e.g., zooplankton) controls? Last but not least, we assessed the 
122temporal changes in the strength of these causal linkages across the studied periods 
123and in how the intensive restoration efforts (IRE) implemented in Lake Taihu after 


1242007 alter their strength. 


125Materials and methods 


126 Data collection 


127 Lake Taihu, the third largest freshwater lake (2,338 km’) in China, is located in a 
128highly developed region characterised by heavy urbanisation and industrialisation in 
129the Yangtze River Delta. The lake has a mean depth of 1.9 m and a catchment area of 
130about 36,500 km’. It is fed by a number of tributaries (>200) that transport massive 
13lamounts of nutrients derived from the watershed into the lake, which is subjected to 
132several rainfall events especially during summer (Yang et al., 2016). Eight sites (i.e., 
133from Meiliang Bay to the centre of the lake, Fig. S1) have been monthly observed by 
134the Taihu Laboratory for Lake Ecosystem Research (TLLER) since 1991, but more 
135regularly since 2000. We compiled monthly observations of Chl (a proxy of 
136phytoplankton biomass), TN, TP, PO4, NO3, NH4 and TIN concentrations, pH, 
137Secchi depth (SD) and zooplankton density (zoo) during the period 2000-2017 from 
138the National Ecosystem Science Data Center (http://www.nesdc.org.cn/). Daily data 
1390n meteorological variables (2000-2017) at the Wuxi Meteorological Station near 
140Lake Taihu, including air temperature (AirTemp), precipitation (Precip) and wind 


141speed (WindSpeed), were collected from the China Meteorological Data Sharing 
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142Service System (http://cdc.cma.gov.cn/) and then calculated monthly average data for 
143use in our analysis. The monthly anomalies (2000-2017) of El Niño and Southern 
144Oscillation (ENSO), Pacific decadal oscillation (PDO) and the Atlantic multidecadal 
145oscillation (AMO) were considered as a proxy of the global climatic oscillation, and 
146oscillation information was obtained from the website of National Weather Service 


147(https://www.weather.gov/). 


148 We applied the monthly averaged data set (2000-2017) on all tested variables, 
149yielding 1728 data records for local environmental variables across 8 sampling sites 
150and 216 data records for the variables of the meteorological and global climatic 
15loscillation (one data record for each month). The limited missing values (< 5%) were 
152interpolated using the “na_seadec” function in the R package of “imputeTS”’ . For 
153each time series, the long-term linear trend was eliminated to ensure stationarity by 
154using the residuals from a linear regression against time. As seasonality could have a 
155confounding influence on the analysis, we deseasonalised the time series by scaling it 
156against the mean and standard deviation of values occurring in the same month 
157according to the equation: (O(¢,)-u:)/o;, where u; is the monthly mean, and ø; is the 
158monthly standard deviation for the i (1, 2, ..., 12) month, and O(t;) is the original time 


159series. 


160Strength of causal linkage and feedback loops 


161 We identified the causal linkages between pairs of time series (i.e., tested 
162variables) using convergent cross mapping (CCM, see Supplementary methods), 


163which is an empirical causality analysis based on Takens’ theorem for dynamical 
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164systems . We quantified the strength of causal linkages as the cross-mapping skill at 
165maximal library size, p(Lma:), 1.e. the greater the p(Lmax), the stronger are the causal 
166effects. For instance, a strong causal effect of PO4 on Chl revealed by CCM indicated 
167that PO4 is a major driver for the dynamics of Chl (magnitude or variability). To 
168eliminate systematic differences (e.g., noise) in cross-mapping skill among the 
169sampling sites, the strength of causal linkages was standardised by dividing linkage 
170strength (LS) by the maximum within each site . The standardised linkage strength 
171(SLS) ranges from 0 to 1, where a value close to | indicates strong causal linkage and 
172vice versa with a value close to 0. We reconstructed causal networks (Fig. 1) and then 
173standardised the strength of causal linkages separately for each site as the distinct sites 


174might not have equivalent dynamics and belong to the same attractor. 


175 The strength of pairwise and triple feedback loops was quantified using Neutel’s 
176loop weight — the geometric mean of SLS for all linkages within a given feedback . In 
177a pairwise feedback (X < Y), the loop weight is the geometric mean of the SLS in 
178both directions (i.e., X — Y and Y — X). Likewise in a triple feedback, the loop 
179weight is the geometric mean of the SLS in a directed cyclic loop (i.e., Type I: X > Y 
180— Z — X; Type I: X — Z — Y — X). To assess the linkage strength of the causal 
18lpathway from global climate oscillation (e.g., ENSO, PDO, AMO) to the local 
182feedback loops via changes in meteorological variables, we calculated the geometric 
183mean of the SLS for all linkages in the pathway (i.e., X — Y — Z). To determine the 


184uncertainty of our estimates in causal strength and loop weight, we calculated their 


20 


185standard errors using the resampling method that reconstructed the sampling 


186distributions from 500 random samples of embedded data points with replacement. 


187Enhancement and regulatory loops for emergent macroecological patterns 


188 We used a generalized least square (GLS) model with a temporal auto- 
189regressive error structure (corArma function) to assess the log-linear relationships 
190between Chl and TN, TP, SD and zooplankton biomass separately for each site. The 
191slope estimated from the GLS model of Chl-TP, Chl-TN, Chl-SD relations and zoo- 
192Chl relations was considered as emerging macroecological patterns of the response of 
193phytoplankton to bottom-up (EMPbottom-up: EMPca-re, EMP cary, EMPecnisp) and top- 
194down (EMPiop-down? EMP o0-cni) controls, respectively. The larger the values of EMP, the 
195stronger the response of phytoplankton to the drivers and vice versa. We tested if there 
196was a significant relationship between the EMP and the strength of causal linkages 
197and feedback loops using Spearman correlation. The feedbacks that have positive 
198relationships with EMP (positive values for EMPcn-rr, EMPcu-rn) and negative 
199relationships (negative values for EMPcn-sp, EMP2oo-cni) could be considered as the 
200enhancement loops of the EMP, promoting the phytoplankton response to the drivers, 


20land vice versa for depressing the phytoplankton response to the drivers. 


202 To further assess the relative importance of these feedbacks for determining the 
203EMP, we used random forests regression (RF) — a bootstrapping classification tree . 
204RF does not only avoid the issue of over-fit (i.e., there are many more responsible 
205variables then for the 8 samples in our study), it also deals efficiently with the spatial 


206autocorrelation among the samples . For each of the four EMPs, we built RF for the 
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207significant causal linkages and feedback loops according to the above Spearman 
208correlation as potential predictors that were randomly resampled to generate 500 un- 
209pruned decision trees. The relative importance of each variable was evaluated 
210according to the increments in mean standard error (MSE) for the RF model 
211predictions. In addition, to assess the relative contribution of enhancement loops and 
212regulatory loops to each EMP, we calculated the percentages of the sum of MSE for 
213all enhancement or regulatory loops dividing them by total MSE for all significant 


214feedback loops with a positive value. 


215 Temporal variations in the strength of causal linkages 


216 We estimated the temporal changes in causal linkages using two complementary 
217approaches. First, we compared the overall strength of causal linkages and feedbacks 
218among the tested variables between the two distinct periods according to the 
219implementation of intensive restoration efforts (IRE) (Pre-IRE: 2000-2007; Post-IRE: 
2202008-2017) by conducting the CCM analysis and measurements of SLS and loop 
221weights separately for each period. Then, we used permutation test (N = 9999) to 
222compare the strength of the two distinct causal linkages or feedback loops across time 
223series (2000-2017) and to test if the same causal linkages or feedback loops differed 


224between the two periods (i.e., Pre-IRE and Post-IRE). 


225 Second, we quantified the time-varying (i.e., dynamic) interaction for each pair 
226of causal linkage (e.g., X — Y). With this, we applied the multivariate S-map method 
227(i.e., a locally weighted multivariate linear regression, see Supplementary methods) to 


228estimate the strength and direction of the interaction between ecological components 
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229by recovering the Jacobian elements at each time point . For each pair of causal 
230linkage, we first calculated the dynamic interaction strength (DIS) as the mean value 
23lof S-map coefficients for each time step across 8 sampling sites and then used 
232Kendall’s z test to assess the significance of temporal trends. The Kendall’s z test was 
233conducted with several time series: 1) whole time series, 2) within seasonal time 
234points (e.g., spring, summer, autumn, winter), 3) within distinct periods (e.g., Pre- 
235IRE: 2000-2007; Post-IRE: 2008-2017). 

236 All statistical tests were performed using R version 4.01 software (R Core Team. 
2372018). CCM and multivariate S-map were performed using the rEDM package 
238 (https://github.com/SugiharaLab/rEDM) in R (v-0.7.5) . The basic ideas of CCM are 
239described in brief animations: tinyurl.com/EDM-intro, and the detailed information to 
240perform EDM analyses can be found at 
24 |https://deepeco.ucsd.edu/resources/#pagecontent. We performed Kendall’s t test using 


242the function mann-kendall in the R package. 


243 Results 


244Causal linkages 


245 Based on CCM analysis (Table S1), we found significant causal linkages with 
246Chl for the local drivers (e.g., SD, pH, NO3, PO4, NH4, TIN, zooplankton)as well as 
247the regional drivers (e.g., AirTemp, Precip, WindSpeed). According to the permutation 
248test for the SLS between the distinct causal linkages, we found that phytoplankton 
249was most sensitive to NH4 and TIN among the tested drivers (Fig. 2A, all P < 0.05) 


250and showed a stronger response to PO4 than to NO3 (Fig. 2A, P < 0.05). Notably, the 
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25laverage strength of the causal effects of local physico-chemical environments (e.g., 
252SD, pH) and zooplankton on Chl was comparable with that of PO4 (Fig. 2A, P > 


2530.05). 


254 In the opposite direction, the SLS of Chl—PO4 was significantly stronger than 
255ChI—NO3 (Fig. 2B, P = 0.027) but not significantly different from Chl—NH4 and 
256ChI-—TIN (Fig. 2B, P > 0.05). The causal effects of phytoplankton on the other local 
257drivers (i.e., SD, pH, zooplankton) were as strong as those on PO4 (Fig. 2B, all P > 
2580.05). However, there was no significant difference between the two directions for the 
259same pairs of variables (Fig. 2A & 2B, e.g., Chl—PO4 vs. PO4—Chl, all P > 0.05) 


260except for a markedly higher NH4—Chl than Chl—NH4 (Fig. 2A & 2B, P < 0.05). 


261 The causal effects of regional drivers on phytoplankton (e.g., AirTemp—Cnhl, 
262Precip—Chl, WindSpeed—Chl) are shown in Fig. 2D, and their effect strengths were 
263significantly lower than those of the local drivers (all P < 0.05). In addition, 
264significant causal effects of global climate oscillation on regional meteorological 
265variables were observed for PDO—WindSpeed, ENSO—AirTemp and 
266ENSO— WindSpeed (Table S3 online), while the strength of their causal pathways to 
267Chl via wind speed was usually stronger than via the other meteorological variables 


268(Fig. 2D, all P< 0.05). 


269 Feedback loops 


270 We quantified the strength of all potential pairwise and triple feedback loops. 
271For the pairwise feedback loops, we found the lowest strength for NO3<>Chl (Fig. 


2722C, all P < 0.05) and the highest for NH4<>Chl (Fig. 2C, all P < 0.05). Notably, the 
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2730ther pairwise feedback loops between phytoplankton and physico-chemical 
274environments (i.e., pH, SD) as well as predators (i.e., zooplankton) were as strong as 
275PO04/TIN/NH4<Chl (Fig. 2C, P > 0.05). Also, the strength of pairwise feedback 
276loops between variables without Chl was assessed (Fig. S3); here TINepH and 
277TIN« zoo had the highest average strength, followed by NH4<>zoo, NO3/NH4<>pH. 
278For the triple feedback loops, there were no significant in-between differences (Fig. 
2792E&2F, all P > 0.05), and the PO4-associated triple feedbacks were as strong as the 


280N-associated triple feedbacks (i.e., NH4, NO3, TIN) ones (Fig. 2E&2F, all P > 0.05) 


281 Enhancement and regulatory loops controlling emergent macroecological patterns 


282 The causal linkages and feedbacks were ubiquitous in Lake Taihu and occurred 
283either as enhancement or regulatory loops, contributing to the emergent 
284macroecological patterns (EMPs) of the phytoplankton response to major drivers (Fig. 
285S4), with 19 significant enhancement loops and 5 regulatory loops for EMP cur (Fig. 
286S4A), 12 enhancement loops and 7 regulatory loops for EMP cut (Fig. S4B), 21 
287enhancement loops and 4 regulatory loops for EMP cusp (Fig. S4C), 5 enhancement 


288loops and 1 regulatory loop for EMP,oo-cni (Fig. S4D). 


289 Random forest analysis revealed that the enhancement loops for EMPs were 
290mainly characterised as local environment- (e.g., nutrients, SD, pH and zooplankton) 
29 associated feedbacks, and the regulatory loops for EMPs were indicative of regional 
292climate- (e.g., wind speed) associated feedbacks (Fig. 3). For these tested feedbacks, a 
293remarkable higher importance of enhancement loops than of regulatory loops was 


294detected in determining EMPen-tp (Fig. 3A), EMPcuisp (Fig. 3C) and EMPcnt-zoo (Fig. 


29 15 
30 


2953D), while the importance of enhancement loops and regulatory loops in affecting 


296EMPcn-r was comparable (Fig. 3B). 


297 Temporal changes in the strength of causal linkages and feedbacks 


298 CCM analysis identified a significant decrease in the strength of 5 causal 
299linkages, 5 causal pathways, 4 pairwise feedbacks and 14 triple feedbacks (Fig. 4). 
300Among these, Windspeed—Chl and AMO-—-Windspeed—Chl were important 
301regulatory loops for EMPbottomup, TIN—Chl and TIN«Chl were important 
302enhancement loops for EMPcu-re, and EMPcn-m, zoo«s>Chl was an important 
303regulatory loop for EMPcnizo. For the triple feedback loops, 
304P04/NH4/TIN—Chl—SD-—PO4/NH4/TIN, NH4/NO3—Chl—zoo—NH4 /NO3 
305were important enhancement loops for EMPcu-re, TIN—>Chl—SD—>TIN was an 
306important enhancement loop for EMPcn-rp and EMPcnisp, PO4—Chl—pH—PO4, 
307NH4/NO3/TIN—Chl—zoo—NH4 /NO3/ TIN were important enhancement loops for 


308EMPcn-sp, SD—pH—Chl—SD was an important enhancement loop for EMPeat-zoo. 


309 S-map analysis revealed that the interactions between the ecological factors are 
310dynamic (i.e., time-varying) rather than static (Fig. 4). We found a significant 
311temporal trend in the dynamic interaction strengths (DIS) of a series of causal 
312linkages (Fig. 4), with decreasing trends for ChlNH4/NO3/P04 and NO3—Chl and 
313increasing trends for ChI—TIN and Precip— Chl. In addition, we observed that DIS 
314differed among seasons (Fig. 6). In spring across the study period (2000-2017), DIS 
315decreased significantly for zoo—Chl. In summer, DIS increased significantly for 


316zoo—Chl and Chl—TIN, while it decreased for PO4—>Chl. In autumn, DIS increased 
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317significantly for zoo/TIN—Cnhl, and decreased for Chl>PO4/NO3/NH4/SD/pH/ and 
318P04/NO3/SD—Chl. In winter, DIS increased significantly for 
319NH4/NO3/AirTemp—Chl and Chl—PO4/NO3/NH4/ pH and decreased for 
320z00/pH/TIN/Precip—Chl. 

321Discussion 

322 We used empirical dynamic modelling (e.g., CCM, S-map) to detect the causal 
323effects and the feedback interactions among the ecological components in Lake Taihu 
324during the last two decades. Our results quantified the casual feedback networks 
325between phytoplankton, water nutrients, physical environments and zooplankton, 
326meteorology as well as global climate oscillation, and quantify the strength of casual 
327linkages and feedback loops. The feedbacks were distinguished into the enhancement 
328and regulatory loops in determining the EMPs of phytoplankton’s response to bottom- 


329up and top-down controls, and their temporal changes were assessed. 


330Causal effects of environmental factors on phytoplankton 


331 We found a greater sensitivity of phytoplankton to NH4 than to NO3 which 
332corresponds well with the fact that phytoplankton generally prefer NH4 due to the 
333lower energetic costs of uptake and assimilation . Previous findings from Lake Taihu 
334also revealed that NH4 was the most important predictor of total phytoplankton 
335biomass and cyanobacteria . The lower causal effects of PO4 than of NH4 on Chl 


336might reflect the lower variation in PO4 (constantly high) in Lake Taihu. 


337 The observed strong effects of zooplankton on Chl suggest that the temporal 


338dynamics of the zooplankton community were encoded in that of phytoplankton, and 
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339we clearly captured the cue of the top-down effects on the phytoplankton in the lake. 
340Consistently, Mao et al. (2020) found that the phytoplankton biomass was positively 
34lassociated with benthivorous and zooplanktivorous fish and negatively with 
342zooplankton biomass, demonstrating that the trophic cascade effects, as fish and 
343zooplankton (43.2%) were comparable with those nutrient effects (40.9%) in 
344controlling the phytoplankton biomass in the lake. In addition, analyses of data from 
345233 Danish lakes, lab and full-scale field experiments highlighted the occurrence of 
346top-down control by fish in shallow lakes . Considering the strong trophic cascade in 
347Lake Taihu, the temporal dynamics of the fish community would be encoded in both 
348the phytoplankton and the zooplankton, and direct and indirect (e.g., trophic cascade) 
349causal effects of the fish community on phytoplankton (top-down forces) would be 


350expected . 


351 Our results elucidated that the local physico-chemical environments (e.g., pH, 
352SD) overall affected Chl to the same extent as PO4 (bottom-up) and zooplankton (top- 
353down), implying that several previous studies on Lake Taihu (focusing on the role of 
354nutrient effects only) might have underestimated the impact of physico-chemical 
355forces on phytoplankton. Besides, the effect strengths of regional drivers (e.g., air 
356temperature, precipitation, wind speed) were significantly lower than those of local 
357drivers (all P < 0.05), which is consistent with previous findings from lab and field 
358studies . Importantly, our results suggest that regional meteorological drivers mediated 
359a cascading effect of the global climate changes on the phytoplankton, as found by 


360other studies . Moreover, except for Chl, significant causal effects of meteorological 
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36lvariables on the other local variables were observed (Fig.S2), supporting the 
362expectation that the drivers from regional to global scale would have a profound 


363influence on the components of the feedback loops . 


364Feedbacks of phytoplankton on environmental factors 


365 Besides the causal effects of environmental drivers on phytoplankton, the 
366comparable feedbacks, both of pairwise ones and triple ones were also documented. 
367The feedbacks were highest for NH4<>Chl, as well as typical bottom-up (nutrients > 
368Chl— zooplankton — nutrients) vs. top-down (nutrients — zooplankton — Chl — 
369nutrients), physico-chemical drivers (nutrients — Chl — pH/SD — nutrients vs. 
370nutrients — pH/SD — Chl — nutrients), characterised physico-chemical 
37lenvironment-mediated predator-prey cycles (pH/SD — Chl — zooplankton — pH/SD 
372vs. pH/SD — zooplankton — Chl — pH/SD). This echo effect supports the idea that 
373phytoplankton has important feedback impacts on the physico-chemical environment, 


374the nutrient recycling, and as well as predators. 


375 The prevalence of the bidirectional causal linkages between pairs of local 
376variables (Fig. S1, Table S2), especially regarding phytoplankton biomass (Chl), 
377provides robust evidence for the hypothesised feedback loops. Phytoplankton first 
378decreases dissolved nutrient abundance by direct uptake, and then return inorganic 
379nutrient back to the water after death; but the rate of nutrient recycling is difficult to 
380measure, and may differ between N and P . The notable causal effects of 
381phytoplankton accompanied with physico-chemical drivers (Fig. S2) on nutrients 


382highlight the importance of water column nutrient recycling processes . Besides, the 


37 19 
38 


383strong feedback strength between zooplankton and especially N (e.g., TIN, NH4) 
384(Fig. S2) emphasizes that the zooplankton community plays an important role in 
385nutrient recycling and phytoplankton growth by excretion . Furthermore, the 
386concurrent changes in the phytoplankton and thus zooplankton community (Ozkan et 
387al, 2014) may entail a highly coupled feedback loops that influence nutrients 


388recycling and stoichiometry of the lake ecosystems. 


389 The feedbacks of Chl to zooplankton (bottom-up) and pH/SD reinforce the 
390comprehensive interactions and thus feedback loops between phytoplankton and the 


391 physico-chemical environment as well as zooplankton (Matsuzaki et al., 2018). In the 


392last four decades, Lake Taihu has experienced major physico-chemical changes 
393including a pronounced rise of pH to 9-10, which may affect both phytoplankton and 
394zooplankton composition and biomass. Algae proliferation fueled by eutrophication 
395could greatly elevate pH and reduce SD, accelerating internal nutrients recycling ). 
396Moreover, the water level has increased 8.7% increase (0.25 m) and the water 
397transparency has decreased 40-60%, which has profound effects on algae blooms and 
398macrophyte loss . Therefore, our results supports that local physico-chemical drivers 
399(e.g., pH, SD, dissolved oxygen, ultraviolet light) as well as predators (i.e., 
400zooplankton) may play a critical mediating role in augmenting or lessening the effects 
401lof phytoplankton on the nutrient recycling . Our results emphasized that these highly 
402interactive components within the constructed causal network create not only efficient 


A403nutrient cycling but also integrative responses to environmental fluctuations. 
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404Enhancement and regulatory loops for emergent macroecological patterns 


405 Causal feedbacks were common in the lake as either enhancement or regulatory 
A06loops contributing to the emergent macroecological patterns (EMPs) of the response 
407of phytoplankton to major drivers. We found that both enhancement and regulatory 
408loops controlled EMPs despite their distinct relative importance. Furthermore, our 
409findings revealed that the algae biomass in the lake might be controlled by both local 
410environment-driven enhancement loops and climate-driven (here wind speed) 
41 lregulatory loops. This suggests that any variation in these local drivers will strengthen 
412(e.g., nutrient inputs, higher pH or turbidity) or weaken (e.g., nutrient reductions, 
413lower pH or increased water transparency) the enhancement loops, which intensified 
414the phytoplankton responses to the bottom-up and top-down controls. In this regard, 
415all loop components act as both “cause” and “effect” variables, and a holistic and 
416integrative view of managing lake ecosystems should be taken. It is noteworthy that 
417the widely used eutrophication mitigation measures in Lake Taihu have greatly 
418reduced the external nutrients input and altered the fish community composition , 
419which (at least for nutrients) weakened the control by the enhancement loops of the 


420phytoplankton biomass, as indicated by our results(Fig. 3&4). 


421 The temporal dynamics of the regional climate strengthen (e.g., increased wind 
422speed) or weaken (e.g., atmosphere stilling) the regulatory loops, thus altering the 
423phytoplankton responses to the tested drivers. Although other variables were not 
424detected as components of the regulatory loops, our results revealed that intrinsically 
425regulatory mechanisms combat the observed self-amplified cycles that sustain algae 


41 21 
42 


426blooms. Moreover, the dominant role of enhancement loops in controlling EMPs 
427(EMPen-te, EMPenisp, EMPcai-zoo) may contribute to the high sensitivity of 
428phytoplankton to variations in TP, water transparency and zooplankton. However, the 
429comparable importance of enhancement loops and regulatory loops in affecting 
430EMPcn-tn may explain the relatively lower sensitivity of phytoplankton to the major 


43 1decreases in TN in the lake. 


432 Temporal changes in the strength of causal linkages and feedbacks 


433 We also quantified the temporal changes in the strength of causal linkages and 
434feedbacks in two dimensions: the overall causal strength based on CCM analysis 
435within the two distinct periods (Pre-IRE: 2000-2007; Post-IRE: 2008-2017) and the 
436dynamic interaction strength (DIS) assessed by the S-map method across the specific 
437time series. Since 2007, a series of intensive restoration efforts (IRE) have been made 
438by the local governments to combat algae blooms and improve water quality . The 
439CCM analysis identified a significant decrease in the strength of several causal 
440feedbacks (Fig. S5, Fig. 6), of which some belong to the enhancement or regulatory 


441loops for the EMPs. 


442 Significantly weaker feedback between phytoplankton and zooplankton were 
443also detected during Post-IRE, which seems reasonable according to long-term 
444observations of trophic interactions in Lake Taihu. Since 2008, biomanipulation by 
445stocking of silver carp and bighead carp to control large-sized phytoplankton (e.g. 
446cyanobacteria) has been conducted, potentially allowing zooplankton to better control 


447the nanophytoplankton . Thus, stronger causal feedback between phytoplankton and 
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448zooplankton would be expected. However, recent findings have indicated that the 
449stocking of phytoplanktivorous fish enhances phytoplankton and cyanobacteria 
450growth and diminishes the zooplankton biomass, the size of cladocerans as well as the 
451zooplankton:phytoplankton ratio , which eventually weakens the grazing pressure of 
452zooplankton on phytoplankton through a trophic cascade. Notably, the S-map results 
453demonstrated a strong seasonal pattern of DIS between zooplankton and 
454phytoplankton with increasing trends during summer and autumn and decreasing 
455trends during winter and spring, which might contribute to a significant seasonality of 


456predator-prey interactions - an area that warrants further studies. 


457 Remarkably, weakened triple feedback loops characterised all the tested local 
458environments (e.g., nutrients, pH, SD, zooplankton), which is greatly associated with 
459the local environment-driven enhancement loops for the EMPs. Therefore, 
460fluctuations in these components will have cascading influences on the efficiency of 
461the enhancement loops in controlling the EMPs even though some of the weakened 
462causal feedbacks (e.g., PO4/pH/SD-associated causal linkages and pairwise 
463feedbacks) are not included as enhancement loops. Overall, this common reduction of 
464the strength of multiple causal interactions and the trophic cascade may give IRE a 


46Scritical role in mitigating the self-amplifying feedback loops supporting algae blooms. 


466 We observed significant decreases in the causal strength (CCM results) of 
467regulatory loops, mainly for climate-associated feedbacks (e.g., wind speed). For 
468Lake Taihu, Deng et al. (2018) revealed that a climatically-modulated decline in wind 


469speed strongly contributed to the variations of Chl, and this decline appears to be 
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470more significant after 2007. In addition, lower wind speed seems to have reduced 
47\non-algal turbidity by nearly 40% and increased the abundance of buoyant 
472cyanobacteria after 2005, providing a “light niche” for phytoplankton and 
A473significantly promoting Chl sensitivity to nutrients . Besides, the significant decrease 
474in the strength of causal pathways (e.g., .ENSO/PDO-Precip—Chl, 
475AMO/PDO—AirTemp—Chl) might also weaken the strength of causal feedback 
476loops through the causal chain from global climate changes to phytoplankton via 
477alterations in precipitation and air temperature. Furthermore, DIS (S-map method) 
478showed increasing trends for AirTemp—Chl during winter, which is consistent with a 
A479previous study identifying a significant association between increased temperature 
480and the sensitivity of Chl to nutrients only in February . These findings indicate that 
48lan alternative mechanism of increased time-lag sensitivity of phytoplankton to 
482nutrients may explain the increasing frequency and magnitude of algae blooms after 
483warmer winters . Overall, the weakened climate (especially wind speed, warming 
484winter)-driven regulatory loops entail increased sensitivity of phytoplankton to local 
A485drivers, which might greatly offset the weaker enhancement loops (external nutrient 


486reduction) by alleviating the phytoplankton response to local drivers . 
487Conclusion 


488 Our results provide a holistic and quantitative evidence for the hypothesised 
489feedback loops and demonstrate how phytoplankton interacts with the tested drivers. 
490The observed common feedback loops reveal that algae proliferation may have a 


491profound impact on the nutrient recycling, the local physico-chemical environment 
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492and zooplankton, while variations in the components of the feedback loop have the 
493potential to lessen or strengthen it. Importantly, we also identified strong cascading 
494causal effects of the global climate change on phytoplankton via changes in regional 
495meteorological drivers (i.e., precipitation and air temperature) that importantly control 
496the dynamics of the feedback loops. Furthermore, our study quantified that the 
497emergent macroecological patterns of the response of phytoplankton to top-down and 
A498bottom-up controls are regulated by two kinds of counterbalance feedback loops: the 
499enhancement loop that strengthens the phytoplankton response and the regulatory 
500loop that suppresses it. Our findings highlight that the constructed causal feedback 
501networks in Lake Taihu loosened significantly after 2007 as a series of local 
502environment-driven enhancement loops underwent a marked decrease in strength 
503during Post-IRE (intense restoration efforts, 2008-2017), indicating that IRE has the 
504potential to combat viscous feedback loops. Moreover, the varying causal effects of 
505global and regional climate changes on the components of the feedback loops will 
506inevitably cause deterioration (i.e., increased sensitivity of phytoplankton to local 
507drivers). Overall, the current IREs initiated to slow down the enhancement loops that 
508trigger over-proliferation of phytoplankton is largely offset by the climate (especially 
509decreasing wind speed, warming winter) and carp stocking, and does thus not suffice 
510to weaken the sensitivity of phytoplankton to the local environment. To counteract the 
511self-amplifying feedback loops, therefore, the present eutrophication mitigation 
512efforts, especially long-term nutrient reduction, should be continued and introduction 


5130f some alternative measures to indirectly regulate the critical components (e.g., pH, 
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514SD, zooplankton) of the loops would also be beneficial. 
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706Fig. 1 The hypothesized feedback loops between phytoplankton biomass (Chl) and 


707local environments (N: nitrogen, P: phosphorus, pH, SD: Secchi depth), zooplankton 
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709climate oscillation (ENSO, PDO, AMO). 

710 
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713 


714Fig. 2 Standardized linkage strength of casual linkages (A and B), loop weight of 
715pairwise feedbacks (C), casual pathway (D) and triangle feedbacks (E and F) between 
716the tested variables across the time series (2000-2017). Chl: chlorophyll a; NO3: 
717nitrate; NH4: ammonium; TIN: total inorganic nitrogen; PO4: orthophosphate; SD: 
718Secchi depth; zoo: zooplankton density. 
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722Fig. 3 Random forest reveals mean predictor importance (percentage of increase of 
723mean square error, Increase in MSE) of each significant casual linkage and feedback 
724in determining the emerging macroecological patterns (EMPs, the estimated slope of 
725the log-log linear relationships) of the phytoplankton response to bottom-up (A: TP, 
726B: TN, C: SD) and top-down (D: zooplankton) controls. The enhancement loops 
727promote phytoplankton response to the drivers with positive relationships with 
728positive EMPs (Chl-TP/TN) and negative relationships with negative EMPs (Chl- 
729SD/zoo), while the regulatory loops depress phytoplankton response to the drivers. 
730The relative contributions of pulled enhancement and regulatory loops were shown. 
731Chl: chlorophyll a; NO3: nitrate, NH4: ammonium; TIN: total inorganic nitrogen; 
732PO4: orthophosphate; SD: Secchi depth; zoo: zooplankton density. 
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734 
735Fig. 4 Significant changes in the standardized linkage strength of causal linkages (A), 


736loop weight of pairwise feedbacks (B) casual pathway (C) and triangle feedbacks (D) 
737between the tested variables during two distinct periods (Pre-IRE: 2000-2007, Post- 
738IRE: 2008-2017). Chl: chlorophyll a; NO3: nitrate; NH4: ammonium; TIN: total 
739inorganic nitrogen; PO4: orthophosphate; SD: Secchi depth; zoo: zooplankton density. 
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743Fig. 5 Time-varying casual interactions with significant temporal trends across time 
744series (monthly from 2000 to 2017) according to Kendall’s t test. Chl: chlorophyll a; 
745NO3: nitrate; NH4: ammonium; PO4: orthophosphate; TIN: total inorganic nitrogen; 
746Precip: precipitation. *: p < 0.05, **: p < 0.01, ***: p < 0.001. 


69 35 
70 


202306.00013v1 


chinaXiv 


747 
748 


z00->Chl AirTemp->Chl pH-=Chl TIN->Chl 


Chi->PO4 Chli->NO3 NH4->Chl 


SD->Chl NO3->Chl Chi->SD Chi->NH4 


ChI->TIN PO4->Chl 


749 


750Fig. 6 Results of Kendall’s t test for the time-varying casual interactions involving 
75\1phytoplankton during four seasons. Chl: chlorophyll a; NO3: nitrate; NH4: 
752ammonium; TIN: total inorganic nitrogen; PO4: orthophosphate; SD: Secchi depth; 
753zoo: zooplankton density, AirTemp: air temperature; Precip: precipitation. *: p < 0.05, 
754**: p < 0.01, ***: p < 0.001, NS: not significant. 
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